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ABSTRACT 

Using precise relative astrometry from the Hubble Space Telescope and the W. M. Keck Telescope, 
we have determined the orbits and masses of the two dynamically interacting satellites of the dwarf 
planet (136108) Haumea, formerly 2003 EL61. The orbital parameters of Hi'iaka, the outer, brighter 
satellite, match well the previously derived orbit. On timescales longer than a few weeks, no Kcplerian 
orbit is sufficient to describe the motion of the inner, fainter satellite Namaka. Using a fully-interacting 
three point-mass model, we have recovered the orbital parameters of both orbits and the mass of 
Haumea and Hi'iaka; Namaka's mass is marginally detected. The data are not sufficient to uniquely 
determine the gravitational quadrupole of the non-spherical primary (described by J2). The nearly 
co-planar nature of the satellites, as well as an inferred density similar to water ice, strengthen the 
hypothesis that Haumea experienced a giant collision billions of years ago. The excited eccentricities 
and mutual inclination point to an intriguing tidal history of significant semi-major axis evolution 
through satellite mean-motion resonances. The orbital solution indicates that Namaka and Haumea 
are currently undergoing mutual events and that the mutual event season will last for the next several 
years. 

Subject headings: comets: general — Kuiper Belt — minor planets — solar system: formation 



1. INTRODUCTION 

The dwarf planet (136108) Haumea, formerly 
2003 EL61, and about 3/4 of other large Kuiper belt 
objects (KBOs) hav e at least one small c lose-in satellite 
(I Weaver et al.l 120061 : IBrown et all 120061 : IBrown fc Suen 
120071 ) . All of these larger KBOs are part of the ex- 
cited Kuiper belt, where the detectable binary frac- 
tion among s maller KBOs is muc h lower, only a 
few percent (jStephens fe Nolll I2006D . In contrast, 
the cold classica l population (inclinati o ns < 5°) ha s 
no large KBOs (|Levison fc Sternl I200H : IBrown! [2008h . 
but prevalent widely separate d binaries with nearly 
equal masses (jNoll et al.l l2008h The differences be- 
tween the types and frequency of Kuiper belt bina- 
ries may point to different binary formation mecha- 
nisms. Small satellites of large KBOs appear to be 
formed by co ll ision, as proposed for the Pluto system 
Canupl I2005t IStern et alj I2006D . Eris and Dysnomia 
Brown fc Schallerl l2007l~ but see iGreenberg fc Barnes! 



2008), and Haume a (Barku me et al.ll2006t IBrown et al.l 
20071 : iFraser fc Brownl I2009D . but smaller KBO binaries 



have more angular momentum than can be generated in 
typical impacts a nd are apparently formed by some other 
mechanism (e.g., [W cidcnschilling 2002; Goldrcich et al.l 
20021 iFunato et al.ll2004t lAstakhov et al.ll2005l : iNesvornvl 
20081) . Both mechanisms of binary formation require 
higher number densities than present in the current 
Kuiper belt, as modeled exp licitly for the Haumea colli- 
sion bv lLevison et alj (|2008f ). 

The collisional origin of Haumea's two satellites — the 
outer, brighter satellite Hi'iaka (SI) and the inner, fainter 
satellite Namaka (S2) — is inferred from several related 
observations. Haumea has a moderate-amplitude light- 
curve and the shortest rotation period (3.9155 hours) 
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among known objects of its size (jRabinowitz et al.ll2006l ). 
The rapid rotation requires a large spin angular mo- 
mentum, as imparted by a large oblique impact. Us- 
ing the mass of Ha umea derived by the orbit of Hi'iaka 
(Bro wn et ail 120051 hereafter B05), assuming Haumea's 
rotation axis is nearly perpendicular to the line-of-sight 
(like the satellites' orbits), and assuming the shape is 
that of a Jacobi ellipsoid (a homogeneous fluid), the pho- 
tometric light curve can be used to d etermine the size, 
shape, albedo, and density of Haume a (IRabinowitz et al.l 
120061 lLacerda fc Jewittll2007l but see lHolsappldl2007f ). It 
is estimated that Haumea is a tri-axial ellipsoid with ap- 
proximate semi-axes of 500 x 750 x 1000 km with a high 
albedo (0.73) and densit y (2.6 g/cm 3 ), as determined by 
iRabinowitz et all (|2006l ). This size and albedo are con - 
sistent with Spitzer radiometry (jStansberrv et all [2008). 
The inferred density is near that of rock and higher than 
all known KBOs implying an atypically small ice frac- 
tion. 

Haumea is also the progenitor of the only known co l- 
lisional family in the Kuiper belt (|Brown et alJ [2007) . 
It seems that the collision that imparted the spin an- 
gular momentum also fragmented and removed the icy 
mantle of the proto-Haumea (thus increasing its den- 
sity) and ejected these fragments into their own he- 
liocentric orbits. The Haumea family members are 
uniquely identified by deep water ice spectra and opti- 
cally neutral color (IBrown et al.ll2007h , flat phase curves 
(IRabinowitz et al.ll2008h. and tight dynamical cluster- 
ing ( Ragozzin e fc Brownl I2007T). The dynamica l clus - 
tering is so significant that iRagozzine fc Brownl (|2007f ) 
were able to correctly predict that 2003 UZ117 and 2005 
CB79 would have deep water ice spe ctra characteristic 
of the Haumea family, as verified by ISchaller fc Brownl 
(2008). The distribution of orbital elements matches the 
unique signature of a collisional family, when resonance 
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diffusion (e.g., iNesvornv fc'R oig 200l|) is taken into ac- 
cou nt. Using this resonance d iffusion as a chronome- 
ter, [Ragozzin e fc Brown! (|2007f ) find that the Haumea 
family-forming collision occurred at least 1 GYr ago and 
is probably primordial. This is consistent with the results 
of iLevison et al.l (|2008[ ). who conclude that the Haumea 
collision is only probable between two scattered-disk ob- 
jects in the early outer solar system when the number 
densities were much higher. 

In this work, we have derived the orbits and masses 
of Haumea, Hi'iaka, and Namaka. In Section 2, we de- 
scribe the observations used to determine precise relative 
astrometry. The orbit-fitting techniques and results are 
given in Section 3. Section 4 discusses the implications 
of the derived orbits on the past and present state of the 
system. We conclude the discussion of this interesting 
system in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

Our data analysis uses observations from various cam- 
eras on the Hubble Space Telescope (HST) and the 
NIRC2 camera with Laser Guide Star Adaptive Objects 
at the W. M. Keck Observatory. These observations are 
processed in different ways; here we describe the general 
technique and below we discuss the individual observa- 
tions. Even on our relatively faint targets (V ~ 21,22), 
these powerful telescopes can achieve relative astrome- 
try with a precision of a few milliarcseconds. The Ju- 
lian Date of observation, the relative astrometric distance 
on-the-sky, and the estimated astrometric errors are re- 
ported in Table [T] 

Observations from Keck are reduced as in B05. Known 
bad pixels were interpolated over and each image di- 
vided by a median flat-field. The images were then 
pair-wise subtracted (from images taken with the same 
filter). The astrometric centroid of each of the visible 
objects is determined by fitting two-dimensional Gaus- 
sians. Converting image distance to on-the-sky astro- 
metric distan ce is achiev e d usi ng the recently derived 
pixel scale of iGhez et al.l {2008), wn0 calibrate the ab- 
solute astrometry of the NIRC2 camera and find a plate 
scale of 0.009963" /pixel (compared to the previously as- 
sumed value of 0.009942" /pixel) and an additional ro- 
tation of 0.13° compared with the rotation info rma- 
tion provided in imag e head ers. IGhez et all (|2008h and 
iHelminiak fc Konackil (|2008l ) find that the plate-scale 
and rotation are stable over the timescale of our observa- 
tions. Error bars are determined from the scatter of the 
measured distances from each individual image; typical 
integration times were about 1 minute. When the inner 
satellite is not detected in individual images, but can be 
seen in the stacked image, then the position is taken from 
the stacked image, after individually rotating, and the er- 
ror bars are simply scaled to the error bars of the outer 
satellite by multiplying by the square root of the ratio 
of signal/noise (^5). The minute warping of the NIRC2 
fields 1 is much smaller than the quoted error bars. 

HST benefits from a known and stable PSF and well- 
calibrated relative astrometry. This allows for precise 
measurements, even when the satellites are quite close 

1 See the NIRC2 Astrometry page at http:// 
www2 . keck.hawaii . edu/inst/nirc2/f orReDoc/ 
post_observing/dewarp/. 



to Haumea. For each of the HST observations, model 
PSFs were generated using Tiny Tim 2 . The model PSFs 
assumed solar colors, as appropriate for Haumea and its 
satellites, and were otherwise processed according to the 
details given in The Tiny Tim User's Guide. All three 
PSFs were then fitted simultaneously to minimize x 2 , 
with errors taken from photon and sky noise added in 
quadrature. Bad pixels and cosmic rays were identi- 
fied by hand and masked out of the \ 2 determinat ion. 
The distortion correction of lAnderson fc Kind (|2003f) for 
WFPC2 is smaller than our error bars for our narrow 
angle astrometry and was not included. Relative on-the- 
sky positions were calculated using the xyad routine of 
the IDL Astro Library, which utilizes astrometry infor- 
mation from the image headers. 

The acquisition and analysis of the satellite images 
taken in 2005 at Keck are described in B05. However, 
there is a sign error in the R. A. Offsets listed in Table 
1 of 2005; the values listed are actually the on-the-sky 
deviations (as visible from their Figure 1). Despite this 
typographical error, the fit of B05 was carried out cor- 
rectly. The observed location s and estima t ed err ors of 
the inner satellite are given in iBrown et all (2006). The 
astrometric positions reported in Table Q] are slightly dif- 
ferent based on a reanalysis of some of the data as well as 
a new plate scale and rotation, discussed above. Based 
on our orbital solution and a reinvestigation of the im- 
ages, we have determined that the May 28, 2005 o bser- 
vation of Namaka reported in IBrown et al.l (|2006f ) was 
spurious; residual long-lived speckles from the adaptive 
optics correction are often difficult to distinguish from 
faint close-in satellites. 

In 2006, HST observed Haumea with the High Resolu- 
tion Camera of the Advanced Camera for Surveys (Pro- 
gram 10545). Two five minute integrations were taken 
at the beginning and end of a single orbit. The raw im- 
ages were used for fitting, requiring distorted PSFs and 
distortion-corrected astrometry. The astrometric accu- 
racy of ACS is estimated to be ~0.1 pixels to which we 
add the photon noise error in the positions of the three 
objects. The high precision of ACS allows for motion to 
be detected between these two exposures, so these errors 
are not based on the scatter of multiple measurements 
as with all the other measurements. 

At the beginning of February 2007, Hubble observed 
Haumea for 5 orbits, obtaining highly accurate positions 
for both satellites (Program 10860). The motion of the 
satellites from orbit to orbit is easily detected, and mo- 
tion during a single orbit can even be significant, so 
we subdivided these images into 10 separate "observa- 
tions". The timing of the observations were chosen to 
have a star in the field of view, from which the Tiny 
Ti m PSF parameter s are m odeled in manner described 
in IBrown fc Truiillol |2004). The observations do not 
track Haumea, but are fixed on the star to get the best 
PSF which is then appropriately smeared for the mo- 
tion of the objects. Even though these observations were 
taken with the Wide Field Planetary Camera — the ACS 
High-Resolution Camera failed only a week earlier — the 
PSF fitting works excellently and provides precise po- 
sitions. Astrometric errors for these observations were 
determined from the observed scatter in positions after 

2 http: / /www. stsci . edu/sof tware/tinytim/t inytim.html. 



Satellites of Haumea 3 



TABLE 1 

Observed Astrometric Positions for the Haumea System 



Julian Date 


Date 


Telescope 


Camera 


Axh 




CTAx H 






A ~ . 
















arcsec 


arcsec 


arcsec 


arcsec 












2453397.162 


2005 


Jan 26 


Keck 


NIRC2 


0.03506 


-0.63055 





.01394 


0.01394 












2453431.009 


2005 


Mar 1 


Keck 


NIRC2 


0.29390 


-1.00626 





.02291 


0.02291 


0.00992 


0.52801 


0.02986 





.02986 


2453433.984 


2005 


Mar 4 


Keck 


NIRC2 


0.33974 


-1.26530 





.01992 


0.01992 












2453518.816 


2005 


May 28 


Keck 


NIRC2 


-0.06226 


0.60575 





.00996 


0.00996 












2453551.810 


2005 


Jun 30 


Keck 


NIRC2 


-0.19727 


0.52106 





.00498 


0.00996 


-0.03988 


-0.65739 


0.03978 





.03978 


2453746.525 


2006 


Jan 11 


HST 


ACS /HRC 


-0.20637 


0.30013 





.00256 


0.00256 


0.04134 


-0.18746 


0.00267 





.00267 


2453746.554 


2006 


Jan 11 


HST 


ACS /HRC 


-0.20832 


0.30582 





.00257 


0.00257 


0.03867 


-0.19174 


0.00267 


0. 


.00267 


2454138.287 


2007 


Feb 6 


HST 


WFPC2 


-0.21088 


0.22019 





.00252 


0.00197 


-0.02627 


-0.57004 


0.00702 





.00351 


2454138.304 


2007 


Feb 6 


HST 


WFPC2 


-0.21132 


0.22145 





.00095 


0.00204 


-0.03107 


-0.56624 


0.00210 





.00782 


2454138.351 


2007 


Feb 6 


HST 


WFPC2 


-0.21515 


0.23185 





.00301 


0.00206 


-0.03009 


-0.55811 


0.00527 





.00564 


2454138.368 


2007 


Feb 6 


HST 


WFPC2 


-0.21402 


0.23314 





.00192 


0.00230 


-0.03133 


-0.56000 


0.00482 





.00663 


2454138.418 


2007 


Feb 6 


HST 


WFPC2 


-0.21705 


0.24202 





.00103 


0.00282 


-0.03134 


-0.54559 


0.00385 





.00376 


2454138.435 


2007 


Feb 6 


HST 


WFPC2 


-0.21449 


0.24450 





.00323 


0.00254 


-0.02791 


-0.54794 


0.00571 


0. 


.00524 


2454138.484 


2007 


Feb 6 


HST 


WFPC2 


-0.21818 


0.25301 





.00153 


0.00224 


-0.02972 


-0.53385 


0.00797 


0. 


.01330 


2454138.501 


2007 


Feb 7 


HST 


WFPC2 


-0.21807 


0.25639 





.00310 


0.00291 


-0.03226 


-0.53727 


0.00531 


0. 


.00400 


2454138.551 


2007 


Feb 7 


HST 


WFPC2 


-0.22173 


0.26308 





.00146 


0.00230 


-0.03429 


-0.53079 


0.00497 





.00582 


2454138.567 


2007 


Feb 7 


HST 


WFPC2 


-0.21978 


0.26791 





.00202 


0.00226 


-0.03576 


-0.52712 


0.00270 


0. 


.00479 


2454469.653 


2008 


Jan 4 


HST 


WFPC2 


0.23786 


-1.27383 





.00404 


0.00824 


-0.02399 


-0.28555 


0.00670 


0. 


.00831 


2454552.897 


2008 


Mar 27 


Keck 


NIRC2 


0.19974 


-0.10941 





.00930 


0.00956 












2454556.929 


2008 


Mar 31 


Keck 


NIRC2 


0.32988 


-0.77111 





.00455 


0.00557 


0.00439 


-0.76848 


0.01239 


0. 


.01280 


2454556.948 


2008 


Mar 31 


Keck 


NIRC2 


0.33367 


-0.77427 





.00890 


0.00753 


0.01363 


-0.76500 


0.01976 





.01252 


2454556.964 


2008 


Mar 31 


Keck 


NIRC2 


0.33267 


-0.77874 





.00676 


0.00485 


0.00576 


-0.77375 


0.01212 





.01283 


2454557.004 


2008 


Mar 31 


Keck 


NIRC2 


0.33543 


-0.78372 





.00404 


0.00592 


0.00854 


-0.77313 


0.01199 





.00897 


2454557.020 


2008 


Mar 31 


Keck 


NIRC2 


0.33491 


-0.78368 





.00374 


0.00473 


0.00075 


-0.76974 


0.00907 


0. 


.01015 


2454557.039 


2008 


Mar 31 


Keck 


NIRC2 


0.33712 


-0.78464 





.00740 


0.00936 


0.00988 


-0.77084 


0.01793 





.01543 


2454557.058 


2008 


Mar 31 


Keck 


NIRC2 


0.33549 


-0.78692 





.00868 


0.00852 


0.01533 


-0.76117 


0.00765 


0. 


.01571 


2454557.074 


2008 


Mar 31 


Keck 


NIRC2 


0.33128 


-0.78867 





.01431 


0.01411 


0.00645 


-0.76297 


0.01639 


0. 


.01390 


2454557.091 


2008 


Mar 31 


Keck 


NIRC2 


0.33687 


-0.79462 





.00803 


0.00717 


0.00708 


-0.76986 


0.01532 


0. 


.00787 


2454593.726 


2008 


May 7 


HST 


NICMOS 


-0.18297 


1.08994 





.00354 


0.00425 


0.00243 


-0.75878 


0.00576 





.00761 


2454600.192 


2008 


May 13 


HST 


WFPC2 


0.10847 


0.17074 





.00508 


0.00427 


-0.02325 


0.19934 


0.00480 





.01161 


2454601.990 


2008 


May 15 


HST 


WFPC2 


0.18374 


-0.13041 





.00729 


0.00504 


-0.02293 


0.50217 


0.00618 


0. 


.00614 


2454603.788 


2008 


May 17 


HST 


WFPC2 


0.24918 


-0.43962 





.00207 


0.00574 


-0.01174 


0.59613 


0.00366 





.00485 


2454605.788 


2008 


May 19 


HST 


WFPC2 


0.29818 


-0.75412 





.00467 


0.00966 


0.00006 


0.29915 


0.00425 


0. 


.00613 



Note. — Summary of observations of the astrometric positions of Hi'iaka (H) and Namaka (N) relative to Haumea. The difference in brightness 
(~6) and orbital planes allow for a unique identi fication of each sate llite without possibility of confusion. The method for obtaining the astrometric 
positions and errors is described in Section[2] and. Brown ct al. (200^). On a few dates, the fainter Namaka was not detected because the observations 
were not of sufficiently deep or Namaka was located within the PSF of Haumea. This data is shown graphically in Figurc[2]and the residuals to the 
fit shown in Figure[3] For reasons described in the text, only the HST data is used to calculate the orbital parameters, which are shown in Table [2] 



subtracting the best fit quadratic trend to the data, so 
that observed orbital motion is not included in the er- 
ror estimate. We note here that combined deep stacks 
of these images revealed no additional outer satellites 
brighter than ~0.25% fractional brightness at distances 
out to about a tenth of the Hill sphere (i.e. about 0.1% of 
the volume where additional satellites would be stable). 

In 2008, we observed Haumea with Keck NIRC2 on 
the nights of March 28 and March 31. The observations 
on March 31 in H band lasted for about 5 hours under 
good conditions, with clear detections of both satellites 
in each image. These were processed as described above. 
Observations where Haumea had a large FWHM were 
removed; about 75% of the data was kept. As with the 
February 2007 HST data, we divided the observations 
into 10 separate epochs and determined errors from scat- 
ter after subtracting a quadratic trend. The motion of 
the outer satellite is easily detected, but the inner satel- 
lite does not move (relative to Haumea) within the errors 
because it is at southern elongation. The March 28 data 
was not nearly as good as the March 31 data due to poor 
weather conditions and only the outer satellite is clearly 
detected. 

In early May 2008, HST observed Haumea using the 
NICMOS camera (Program 11169). These observations 
were processed as described above, though a few images 



with obvious astrometric errors (due to the cosmic rays 
which riddle these images) were discarded. These are th e 
same observations discussed bv lFraser &: Brownl {2009). 

In mid-May 2008, we observed Haumea at five epochs 
using the Wide Field Planetary Camera (WFPC2), over 
the course of 8 days (Program 11518). Each of these 
visits consisted of four ^10 minute exposures. These 
data, along with an observation in January 2008, were 
processed as described above. Although we expect that 
some of these cases may have marginally detected motion 
of the satellites between the four exposures, ignoring the 
motion only has the effect of slightly inflating the error 
bars for these observations. Namaka was too close to 
Haumea (< 0.1") to observe in the May 12, 2008 image, 
which is not used. 

The derived on-the-sky relative astrometry for each 
satellite, along with the average Julian Date of the ob- 
servation and other information are summarized in Ta- 
ble [TJ These are the astrometric data used for orbit 
fitting in this paper. In earlier attempts to determine 
the orbit of Namaka, we also obtained other observa- 
tions. On the nights of April 20 and 21, 2006, we ob- 
served Haumea with the OSIRIS camera and LGSAO at 
Keck. Although OSIRIS is an integral-field spectrome- 
ter, our observations were taken in photometric mode. In 
co-added images, both satellites were detected on both 
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nights. We also received queue-scheduled observations 
of Haumea with the NIRI camera on Gemini and the 
LGSAO system Altair. In 2007, our Gemini program re- 
sulted in four good nights of data on April 9 and 13, May 
4, and June 5. In 2008, good observations were taken on 
April 20, May 27, and May 28. In each of the Gemini 
images, the brighter satellite is readily found, but the 
fainter satellite is often undetectable. 

The accuracy of the plate scale and rotation required 
for including OSIRIS and Gemini observations is un- 
known, so these data are not used for orbit determina- 
tion. We have, however, projected the orbits derived 
below to the positions of all known observations. The 
scatter in the Monte Carlo orbital suites (described be- 
low) at the times of these observations is small compared 
to the astrometric error bars of each observation, imply- 
ing that these observations are not important for improv- 
ing the fit. Predicted locations do not differ significantly 
from the observed locations, for any obse rvation of which 
we are awa re, includ i ng th ose reported in Barkume et aTl 
(2006) and lLacerdal (|2008f ). 

Using these observations, we can also do basic relative 
photometry of the satellites. The brightness of the satel- 
lites was computed from the height of the best-fit PSFs 
found to match the May 15, 2008 HST WFPC2 observa- 
tion. Based on the we ll-known period and phase of the 
light curve of Haumea (jLacerda et alJl2008t D. Fabry cky, 
pers. comm.), Haumea was at its faintest during these 
observations and doesn't change significantly in bright- 
ness. Hi'iaka was found to be ~10 times fainter than 
Haumea and Namaka ~3.7 times fainter than Hi'iaka. 

3. ORBIT FITTING AND RESULTS 

The orbit of Hi'iaka and mass of Haumea were origi- 
nally determined by B05. From three detections of Na- 
maka, iBrown et al.l (|2006f) estimated three possible or- 
bital periods around 18, 19, and 35 days. The ambiguity 
resulted from an under-constrained problem: at least 4-5 
astrometric observations are required to fully constrain 
a Keplerian orbit. Even after additional astrometry was 
obtained, however, no Keplerian orbit resulted in a rea- 
sonable fit, where, as usual, goodness-of-fit is measured 
by the \ 2 statistic, and a reduced \ 2 01 order unity is 
required to accept the orbit model. By forward inte- 
gration of potential Namaka orbits, we confirmed that 
non-Keplerian perturbations due to Hi'iaka (assuming 
any reasonable mass) causes observationally significant 
deviations in the position of Namaka on timescales much 
longer than a month. Therefore, we expanded our orbital 
model to include fully self-consistent three-body pertur- 
bations. 

3.1. Three Point-Mass Model 

Determining the orbits and masses of the full system 
requires a 15-dimensional, highly non-linear, global % 2 
minimization. We found this to be impractical without 
a good initial guess for the orbit of Namaka to reduce the 
otherwise enormous parameter space, motivating the ac- 
quisition of multiple observations within a short enough 
timescale that Namaka's orbit is essentially Keplerian. 
Fitting the May 2008 HST data with a Keplerian model 
produced the initial guess necessary for the global mini- 
mization of the fully-interacting three point-mass model. 
The three point-mass model uses 15 parameters: the 



masses of the Haumea, Hi'iaka, and Namaka, and, for 
both orbits, the osculating semi-major axis, eccentricity, 
inclination, longitude of the ascending node, argument 
of periapse, and mean anomaly at epoch HJD 2454615.0 
(= May 28.5, 2008). All angles are defined in the J2000 
ecliptic coordinate system. Using these orbital elements, 
we constructed the Cartesian locations and velocities at 
this epoch as an initial condition for the three-body in- 
tegration. Using a sufficiently small timestep (~300 sec- 
onds for the final iteration), a FORTRAN 90 program 
integrates the system to calculate the positions relative 
to the primary and the positions at the exact times of 
observation are determined by interpolation. (Observa- 
tion times were converted to Heliocentric Julian Dates, 
the date in the reference frame of the Sun, to account 
for light-travel time effects due to the motion of the 
Earth and Haumea, although ignoring this conversion 
does not have a significant effect on the solution.) Using 
the JPL HORIZONS ephemeris for the geocentric posi- 
tion of Haumea, we vectorially add the primary-centered 
positions of the satellites and calculate the relative astro- 
metric on-the-sky positions of both satellites. This model 
orbit is then compared to the data by computing \ 2 i n 
the normal fashion. We note here that this model does 
not include gravitational perturbations from the Sun or 
center-of-light/center-of-mass corrections, which are dis- 
cussed below. 

Like many multi-dimensional non-linear minimization 
problems, searching for the best-fitting parameters re- 
quired a global minimization algorithm to escape the 
ubiquitous local minima. Our algorithm for finding the 
global minimum starts with thousands of local minimiza- 
tions, executed with the Levenberg-Marquardt algorithm 
mpf it 3 using numerically-determined derivatives. These 
local minimizations are given initial guesses that cover 
a very wide range of parameter space. Combining all 
the results of these local fits, the resultant parameter vs. 
X 2 plots showed the expected parabolic shape (on scales 
comparable to the error bars) and these were extrapo- 
lated to the their minima. This process was iterated 
until a global minimum is found; at every step, random 
deviations of the parameters were added to the best-fit 
solutions, to ensure a full exploration of parameter space. 
Because many parameters are highly correlated, the abil- 
ity to find the best solutions was increased significantly 
by adding correlated random deviations to the fit pa- 
rameters as determined from the covariance matrix of 
the best known solutions. We also found it necessary 
to optimize the speed of the evaluating of \ 2 from the 
15 system parameters; on a typical fast processor this 
would take a few hundredths of a second and a full local 
minimization would take several seconds. 

To determine the error bars on the fit parameters, 
we use a Monte C arlo technique (B05), as suggested in 
iPress et all ()1992l ). Synthetic data sets are constructed 
by adding independent Gaussian errors to the observed 
data. The synthetic data sets are then fit using our global 
minimization routine, resulting in 86 Monte Carlo real- 
izations; four of the synthetic data-sets did not reach 
global minima and were discarded, having no significant 
affect on the error estimates. One-sigma one-dimensional 

An IDL routine available at http:// 
www. physics .wise . edu/ ^craigm/idl/f itting.html. 
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Fig. 1. — Relative positions of the satellites as viewed from Earth. 
The outer orbit corresponds to the brighter Hi'iaka and the inner 
orbit corresponds to the fainter Namaka. In the center is Haumea, 
drawn to scale, assuming an ellipsoid cross-section of 500 x 1000 km 
IIRabino witz ct al. 2006) with the long axis oriented North-South. 
The apparent orbit changes due to parallax and three-body effects; 
this is the view near March 2008. See Figure[2]for model and data 
positions throughout the observation period (2005-2008). 



error bars for each parameter are given by the standard 
deviation of global-best parameter fits from these syn- 
thetic datasets. For each parameter individually, the dis- 
tributions were nearly Gaussian and were centered very 
nearly on the best-fit parameters determined from the 
actual data. The error bars were comparable to error 
bars estimated by calculating where y 2 increas ed by 1 
from the global minimum fsee lPress et a,l.|[l992f ) . 

First, we consider a solution using only the observa- 
tions from HST. Even though these are taken with differ- 
ent instruments (ACS, NICMOS, and mostly WFPC2), 
the extensive calibration of these cameras allows the di- 
rect combination of astrometry into a single dataset. 
The best-fit parameters and errors are shown in Table 
[SJ The reduced chi-square for this model is x 2 e d = 0-64 
(x 2 = 36.4 with 57 degrees of freedom). The data are 
very well-fit by the three point mass model, as shown in 
Figures [5] and [31 A reduced x 2 l ess than 1 is an indi- 
cation that error bars are overestimated, assuming that 
they are independent; we note that using 10 separate 
"observations" for the Feb 2007 data implies that our 
observations are not completely independent. Even so, 
X 2 ed values lower than 1 are typical for this kin d of as- 
trometric orbit fitting (e.g.. iGrundv et aT1l2008D . Each 
of the fit parameters is recovered, though the mass of 
Namaka is only detected with a 1.2-cr significance. Na- 
maka's mass is the hardest parameter to determine since 
it requires detecting minute non-Keplerian perturbations 
to orbit of the more massive Hi'iaka. The implications 
of the orbital state of the Haumea system are described 
in the next section. We also list in Table [3J the initial 
condition of the three-body integration for this solution. 

The HST data are sufficient to obtain a solution for 
Hi'iaka's orbit that is essentially the same as the orbit 
obtained from the initial Keck data in B05. Neverthe- 
less, the amount and baseline of Keck NIRC2 data is 



useful enough to justify adding this dataset to the fit. 
Simply combining these datasets and searching for the 
global minimum results in a significant degradation in 
the fit, going from a reduced x 2 of 0.64 to a reduced \ 2 of 
~1.10, although we note that this is still an adequate fit. 
Adding the Keck data has the effect of generally lowering 
the error bars and subtly changing some of the retrieved 
parameters. Almost all of these changes are within the 
~l-er error bars of the HST only solution, except for the 
mass estimate of Namaka. Adding the Keck data results 
in a best-fit Namaka mass a factor of 10 lower than the 
HST data alone. The largest mass retrieved from the 
entire Monte Carlo suite of solutions to the HST+Keck 
dataset is ~ 8 x 10 17 kg, i.e. a Namaka/Haumea mass 
ratio of 2 x 10~ 4 , which is inconsistent with the bright- 
ness ratio of ^0.02, for albedos less than 1 and densi- 
ties greater than 0.3 g/cc. However, this solution as- 
sumes that the Ke ck NIRC2 absolute astrometry (based 
on the solution of Gh ez et all (|2008l ). which is not di- 
rectly cross-calibrated with HST) is perfectly consistent 
with HST astrometry. In reality, a small difference in 
the relative plate scale and rotation between these two 
telescopes could introduce systematic errors. Adding fit- 
ted parameters that adjust the plate scale and rotation 
angle does not help, since this results in over-fitting, as 
verified by trial fitting of synthetic datasets. We adopt 
the HST-only solution, keeping in mind that the nominal 
mass of Namaka may be somewhat overestimated. 

Using the Monte Carlo suite of HST-only solutions, 
we can also calculate derived parameters and their er- 
rors. Using Kepler's Law (and ignoring the other satel- 
lite), the periods of Hi'iaka and Namaka are 49.462 ± 
0.083 days and 18.2783 ± 0.0076 days, respectively, with 
a ratio of 2.7060 ± 0.0037, near the 8:3 resonance. The 
actual mean motions (and resonance occupation) will be 
affected by the presence of the other satellite and the 
non-spherical nature of the primary (discussed below) . 

The mass ratios of the satellite to Haumea are 0.00451 
± 0.00030 and 0.00051 ± 0.00036, respectively and the 
Namaka/Hi'iaka mass ratio is 0.116 ± 0.086. The mu- 
tual inclination of the two orbits is <f> = 13.41° ± 0.08°, 
where the mutual inclination is the actual angle be- 
tween the two orbits, given by cos0 = cosi# cosijv + 
smitf sinijv cos(fi# — f2jy), where i and Q are the in- 
clination and longitude of ascending node. The origin of 
this significantly non-zero mutual inclination is discussed 
in Section 14.3.21 The mean longitude, A = £1 + oj + M, 
is the angle between the reference line (J2000 eclip- 
tic first point of Ares) and is determined well; the er- 
rors in the argument of periapse (u>) and mean anomaly 
(M) shown in Tabel [2] are highly anti-correlated. Our 
Monte Carlo results give \h — 153.80 ±0.34 degrees and 
\n = 202.57 ± 0.73 degrees. Finally, under the nom- 
inal point-mass model, Namaka's argument of periapse 
changes by about -6.5° per year during the course of the 
observations, implying a precession period of about 55 
years; the non-Keplerian nature of Namaka's orbit is de- 
tected with very high confidence. 

3.2. Including the J2 of Haumea 

The non-spherical nature of Haumea can introduce ad- 
ditional, potentially observable, non-Keplerian effects. 
The largest of these effects is due to the lowest-order 
gravitational moment, the quadrupole term (the dipole 
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Fig. 2. — Observed positions and model positions of Hi'iaka and Namaka. From top to bottom, the curves represent the model on-the-sky 
position of Hi'iaka in the a;-direction (i.e. the negative offset in Right Ascension), Hi'iaka in the ^/-direction (i.e. the offset in Declination), 
Namaka in the ^-direction, and Namaka in the ^/-direction, all in arcseconds. Points represent astrometric observations as reported in Table 
[T] Error bars are also shown as gray lines, but are usually much smaller than the points. The three-point mass model shown here is fit to 
the HST-data only, with a reduced Xred °^ 0-64. The residuals for this solution are shown in Figure [3] Note that each curve has its own 
scale bar and that the curves are offset for clarity. The model is shown for f260 days, starting on HJD 2453297.0, ~f00 days before the 
first observation and ending just after the last observation. Visible are the orbital variations (~49.5 days for Hi'iaka and ~f8.6 days for 
Namaka), the annual variations due to Earth's parallax, and an overall trend due to a combination of Haumea's orbital motion and the 
precession of Namaka's orbit. 



moment is in the center of mass fra me), described by 
J 2 (see, e.g.. lMurrav h Dermottl 2000). Haumea rotates 
over 100 times during a single orbit of Namaka, which or- 
bits quite far away at ~35 primary radii. To lowest order, 
therefore, it is appropriate to treat Haumea as having 
an "effective" time-averaged J 2 . Using a code provided 
by E. Fahnestock, we integrated trajectories similar to 
Namaka's orbit around a homogeneous rotating tri-axial 
ellipsoid and have confirmed that the effective J 2 model 
deviates from the full model by less than half a milliarc- 
second over three years. 

The value of the effective J 2 (= — C20) for a rotating 
homog eneous tri-axial ellipsoid was derived by IScheered 
(fl99l : 

J 2 R 2 = ^{a 2 +(3 2 - 2 7 2 ) ~ 1.04 x 10 n m 2 (1) 

where a, (3, and 7 are the tri-axial radii and the numerical 
value corresponds to a (49 8 x 759 x 980) km ellips oid as 
inferred from photometry ([Rabinowitz et ail l2006). Wc 



note that the physical quantity actually used to deter- 
mine the orbital evolution is J 2 R 2 ; in a highly triaxial 
body like Haumea, it is not clear how to define R, so 
using J 2 R 2 reduces confusion. If R is taken to be the 
volumetric effective radius, then R ~ 652 km and the 
J 2 ~ 0.244. Note that the calculation and use of J 2 
implicitly requires a definition of the rotation axis, pre- 
sumed to be aligned with the shortest axis of the ellip- 
soid. 

Preliminary investigations showed that using this value 
of J2-R 2 implied a non-Keplerian effect on Namaka's orbit 
that was smaller, but similar to, the effect of the outer 
satellite. The primary observable effect of both J 2 and 
Hi'iaka is the precession of a pses and nodes of Namaka 's 
eccentric and inclined orbit ([Murray fc DermottJ l2000n . 
When adding the three relevant parameters — J2-R 2 and 
the direction of the rotational axis 4 — to our fitting pro- 

4 When adding J2, our three-body integration was carried out 
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TABLE 2 

Fitted Parameters of the Haumea System 



Object 


Parameter 


Value 




Error 


Units 


Haumea 


Mass 


4.006 


± 


0.040 


10 kg 


Hi'iaka 


Mass 


1.79 


± 


0.11 


10 kg 




Semi-major axis 


49880 


± 


198 


km 




Eccentricity 


U.Uolo 


± 


n novo 






Inclination 


126.356 


± 


0.064 


degrees 




Longitude of ascending node 


206.766 


± 


0.033 


degrees 




Argument of periapse 


154.1 


± 


5.8 


degrees 




Mean anomaly 


152.8 


± 


6.1 


degrees 


Namaka 


Mass 


1.79 


± 


1.48 


10 18 kg 




Semi-major axis 


25657 


± 


91 


km 




Eccentricity 


0.249 


± 


0.015 






Inclination 


113.013 


± 


0.075 


degrees 




Longitude of ascending node 


205.016 


± 


0.228 


degrees 




Argument of periapse 


178.9 


± 


2.3 


degrees 




Mean anomaly 


178.5 


± 


1.7 


degrees 



Note. — Orbital parameters at epoch HJD 2454615.0. The nominal values 
are from the best fit to the HST data, while the (often correlated) error bars 
are the standard deviation of fitted values returned from a Monte Carlo suite 
of 86 datasets. These are the osculating orbital elements at this epoch; due 
to the three-body interactions, these values (especially the angles) change over 
the timcscalc of observations. All angles are referenced to the J2000 ecliptic 
coordinate system. See Table [3] for the Cartesian positions of the two satellites 
at this epoch. 



TABLE 3 

State Vector for the Haumea System 



Object x (m) y (m) z (m) v x (m/s) v v (m/s) v z (m/s) 

Hi'iaka -18879430 -36260639 -32433454 60.57621 1.85403 -34.81242 
Namaka -28830795 -13957217 -1073907 16.07022 -26.60831 -72.76764 



Note. — Cartesian position and velocity of Haumea's satellites relative to Haumca 
in the J2000 ecliptic coordinate system at epoch HJD 2454615.0 corresponding to the 
best-fit orbital parameters shown in Table [2] 



cedure, we found a direct anti-correlation between J2R 2 
and the mass of Hi'iaka, indicating that these two param- 
eters are degenerate in the current set of observations. 
The value of the reduced % 2 was lowered significantly by 
the addition of these parameters, with the F-test return- 
ing high statistical significance. However, the best fits 
placed the satellites on polar orbits. We have verified 
that the fitted values of J2R 2 and the spin pole perpen- 
dicular to the orbital poles is due to over-fitting of the 
data. We generated simulated observations with the ex- 
pected value of J2R 2 and with the satellites in nearly 
equatorial orbits. Allowing the global fitter to vary all 
the parameters resulted in an over-fitted solution that 
placed the satellites on polar orbits. Hence, allowing J2 
and the spin pole to vary in the fit is not justified; the 
effect of these parameters on the solution are too small 
and/or too degenerate to detect reliably. Furthermore, 
since the model without these parameters already had a 
reduced \ 2 l ess than 1, these additional parameters were 
not warranted in the first place. 

It is interesting, however, to consider how includ- 
ing this effect would change the determination of the 
other parameters. We therefore ran an additional set 
of models with a fixed J2R 2 and fixing the spin pole 
(more accurately, the axis by which J2 is defined) as the 

in the frame of the primary spin axis and then converted back to 
ecliptic coordinates. 



mass- weighted orbital pole of Hi'iaka and Namaka; since 
Hi'iaka is ~10 times more massive than Namaka, this 
puts Hi'iaka on a nearly equatorial orbit (i ~ 1°), as 
would be expected from collisional formation. Holding 
J 2 R 2 fixed at 1.04 x 10 11 m 2 , we reanalyzed the HST 
dataset using our global fitting routine. As expected, 
none of the parameters change by more than 1-er, except 
for the mass of Hi'iaka, which was reduced by almost 
30% to - 1.35 x 10 19 kg. (In fits where J 2 R 2 was al- 
lowed to vary, the tradeoff between J2R 2 and Hi'iaka's 
mass was roughly linear, as would be expected if the sum 
of these effects were forced to match the observationally 
determined precession of Namaka.) The data were well 
fit by the non-point mass model, with the forced J2R 2 
and spin pole solution reaching a global reduced x 2 of 
0.72. Since Haumea's high- amplitude light-curve indi- 
cates a primary with a large quadrupole component, the 
nominal mass of Hi'iaka in the point-mass case is almost 
certainly an overestimate of its true mass. More data 
will be required to disentangle the degeneracy between 
the mass of Hi'iaka and the J2 of Haumca. Including J2 
and/or the Keck data do not improve the estimates of 
Namaka's mass. 

4. IMPLICATIONS OF ORBITAL SOLUTIONS 

Taking the orbital solutions derived in the previous 
section, we can begin to answer questions relevant to the 
formation and evolution of this unique satellite system. 
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Fig. 3. — Normalized residuals of the three-point mass fit to 
HST-data only. Plotted is (Ax mod — Ax ^ a )/<j& x versus (Ay mod — 
Aj/ t, s )/o"Ay f° r Hi'iaka (diamonds) and Namaka (circles). Points 
that lie within the circle indicate where the model and observations 
vary by less than 1 error bar. See also Figure [2] The residuals 
are roughly evenly spaced and favor neither Hi'iaka nor Namaka, 
implying that there are no major systematic effects plaguing the 
three-body fit. As reported in the text this solution has a reduced 
X 2 Icd of 0.64. 



4.1. Mutual Events and Satellite Sizes 

According to the orbit solution, the Haumea system 
is currently undergoin g mutual events, as reported in 
(jFabrvckv et all 120081 ). (Th is is also true using the 
other orbit solutions, e.g. HST+Keck, with or without 
J2O Using the known orbit, the angle between Namaka, 
Haumea, and the Earth (in the case of occultations) or 
the Sun (in the case of shadowing) falls well below the 
~13 milliarcseconds (~500 km) of the projected shortest 
axis of Haumea. Observing multiple mutual events can 
yield accurate and useful measurements of several system 
properties as shown by the r esults of the Pluto-Charo n 
mutual event season (e.g. iBinzel fc Hubbarl fl997h . 
The depth of an event where Namaka occults Haumea 
leads to the ratio of albedos and, potentially, a surface 
albedo map of Haumea, which is known to exhibit color 
variations as a functi on of rotational phase, indicative o f 
a variegated surface (|Lacerda et al.ll200ll lLace rda 2008) . 
Over the course of a single season, Namaka will traverse 
several chords across Haumea allowing for a highly 
accurate measureme nt of Haumea's size, sh ape, and spin 
pole direction (e.g.. iDescamps et al.l 2008). The precise 
timing of mutual events will also serve as extremely 
accurate astrometry, allowing for an orbital solution 
much more precise than reported here. We believe 
that incorporating these events into our astromctric 
model will be sufficient to independently determine the 
masses of all three bodies and J2R 2 • Our solution also 
predicts a satellite-satellite mutual event in July 2009 — 
the last such event until the next mutual event season 
begins around the year 2100. Our knowledge of the 
state of the Haumea system will improve significantly 
with the observation and analysis of these events. See 
http : //web . gps . caltech. edu/^mbrown/2003EL61/mutual 



for up-to-date information on the Haumea mutual events. 
Note that both the mutual events and the three-body 
nature of the system are valuable for independently 
checking the astrometric analysis, e.g. by refining plate 
scales and rotations. 

Using the best-fit mass ratio and the photometry de- 
rived in Section [21 we can estimate the range of albedos 
and densities for the two satellites. The results of this 
calculation are shown in Figured] The mass and bright- 
ness ratios clearly show that the satellites must either 
have higher albedos or lower densities than Haumea; the 
difference is probably even more significant than shown 
in Figure |H since the nominal masses of Hi'iaka and Na- 
maka are proba bly overestimated (see Section [3| . The 
similar spectral (jBarkume et alJfeOOrjf ) and photometric 
(jFraser fc Brownl l2009f ) properties of Haumea, Hi'iaka, 
and Namaka indicate that their albedos should be sim- 
ilar. Similar surfaces are also expected fr om rough ca l- 
culatio ns of ejecta exchange discussed by ISternl (2009), 
though iBenecchi et alJ (|2008h provide a contrary view- 
point. If the albedos are comparable, the satellite den- 
sities indicate a mostly water ice composition (p w 1.0 
g/cc). This lends support to the hypothesis that the 
satellites are formed from a collisional debris disk com- 
posed primarily of water ice from the shattered mantle 
of Haumea. This can be confirmed in the future with a 
direct measurement of Namaka's size from mutual event 
photometry. Assuming a density of water ice, the es- 
timated radii of Hi'iaka and Namaka are ^160 km and 
~80 km, respectively. 

4.2. Long-term Orbital Integrations 

It is surprising to find the orbits in an excited state, 
both with non-zero eccentricities and with a rather large 
mutual inclination. In contrast the regular satellite 
systems of the gas gian ts, the satellites o f Mars, the 
three satellites of Pluto (jTholen et al.ll2008|). and aster- 
oid tr iple systems with well known orbits (jMarchis et alJ 
120051) are all in nearly circular and co-planar orbits. In 
systems of more than one satellite, perturbations be- 
tween the satellites produce forced eccentricities and in- 
clinations that will remain even with significant damping. 
If the excited state of the Haumea system is just a reflec- 
tion of normal interactions, then there will be small free 
eccentricities and inclinations, which can be estimated by 
integrating the system for much longer than the preces- 
sion timescales and computing the time average of these 
elements. Using this technique, and exploring the entire 
Monte Carlo suite of orbital solutions, we find that the 
free eccentricity of Hi'iaka is ~0.07, the free eccentric- 
ity of Namaka is ~ 0.21, and the time-averaged mutual 
inclination is ^12.5°. Non-zero free eccentricities and in- 
clinations imply that the excited state of the system is 
not due to satellite-satellite perturbations. These inte- 
ntions were calculated u sing the n-body code SyMBA 
Levison &: Duncan! Il994h using the regularized mixed 
variable symplectic meth od based on the mapping by 
IWisdom fc Holnianl (|1991f ). Integration of all the Monte 
Carlo orbits showed for ~2000-years showed no signs of 
instability, though we do note that the system chaotically 
enters and exits the 8:3 resonance. The orbital solutions 
including J2R 2 — 1.04 x 10 11 m 2 were generally more 
chaotic, but were otherwise similar to the point-mass in- 
tegrations. 




Satelli tes of Haumca 

i4oi"i""i""i""i""i! — 




0.5 1.0 1.5 2.0 2.5 3.0 
Density (g/cc) 



0.5 1.0 1.5 2.0 2.5 3.0 
Density (g/cc) 



Fig. 4. — Relationship between radius, density, and albedo for 
Hi'iaka (left) and Namaka (right). A range of possible albedos 
and densities can reproduce the determined mass and brightness 
ratios of Hi'iaka and Namaka, which are assumed to be spheri- 
cal. The solid lines show the relationship for the nominal masses, 
reported in Table [2] with dotted lines showing the l-a mass er- 
ror bars. Note that the mass of both Hi'iaka and Namaka may 
be overestimated (due to insufficient data, see S ection |3j. The 
albedo and density of the Rabinowitz ct al. (2006) edge-on model 
for Haumea are shown by dotted lines. Both satellites must have 
lower de nsities and/or higher albedos than Hau mea. The similar 
spect ral IIBarkume et al.ll2006l) and photometric (Fraser fc Brownl 
2009) properties of Haumea, Hi'iaka, and Namaka indicate that 
their albedos should be similar. Under the assumption that the 
satellites have similar albedos to Haumea, the densities of the 
satellites indicate that they are primarily composed of water ice 
(p ~ 1.0 g/cc). Low satellite densities would bolster the hypothe- 
sis that the satellites formed from the collisional remnants of the 
water ice mantle of the differentiated proto-Haumea. Observation 
of Haumea-Namaka mutual events will allow for much more precise 
and model-independent measurements of Namaka's radius, density, 
and albedo. 



These integrations did not include the effect of the 
Sun, which adds an additional minor torque to the sys- 
tem that is negligible (Af2 ~ 10 -5 degrees) over the 
timescale of observations. The effects of the sun on the 
satellite orbits on long-time scales were not investigated. 
While the relative inclination between the satellite or- 
bits and Haumca's heliocentric orbit (~119° for Hi'iaka 
and ~1 05° for Namaka) places this system i n the regime 
where the Kozai ef fect can be important (jKozail fl962l : 
iPerets fc Naodl2008[ ) , the interactions between the satel- 
lites are strong enough to suppress weak Kozai oscilla- 
tions due to the Sun, w hich are only active in the abs ence 
of other perturbations (jFabrvckv fc Trem ainc 20Q3). 

We did not include any correction to our solution for 
possible differences between the center of light (more pre- 
cisely, the center of fitted PSF) and cen ter of mass of 
Haum ea. This may be important since lLacerda et al.l 
( 2008) find that Haumea's two-peaked light curve can be 
explained by a dark red albedo feature, which could po- 
tentially introduce an systematic astrometric error. The 
February 2007 Hubble data and March 2008 Keck data, 
both of which span a full rotation, do not require a 
center-of-light/center-of-mass correction for a good fit, 
implying that the correction should be smaller than ~2 
milliarcseconds (i.e. ~70 km). Examination of the all the 
astrometric residuals and a low reduced-;^ 2 confirm that 
center-of-light/center-of-mass corrections are not signif- 
icant at our level of accuracy. For Pluto and Charon, 



albedo features can result in spurious orbital astrometry 
because Pluto and Charon are spin-locked; this is not the 
case for Haumea. For Keck observations where Namaka 
is not detected (see Table [T]), it is usually because of low 
signal-to-noise and Namaka's calculated position is not 
near Haumea. In the cases where Namaka's light con- 
taminates Haumea, the induced photocenter error would 
be less than the observed astrometric error. 

4.3. Tidal Evolution 

All of the available evidence points to a scenario for the 
formation of Haumea's satellites similar to the formation 
of the Earth's moon: a large oblique collision created a 
disk of debris composed mostly of the water ice man- 
tle of a presumably differentiated proto-Haumea. Two 
relatively massive moons coalesced from the predomi- 
nantly water-ice disk near the Roche lobe. Interestingly, 
in studying the formation of Earth's Moon , about one 
third of the simulations of llda et al.l (|1997t ) predict the 
formation of two moonlets with the outer moonlet ~10 
times more massive the inn er moonlet. Alt hough the 
disk accretion model used bv llda et al.1 (|1997h made the 
untenable assumption that the remnant disk would im- 
mediately coagulate into solid particles, the general idea 
that large disks with sufficient angular momentum could 
result in two separate moons seems reasonable. Such 
collisional satellites coagulate near the Roche lobe (e.g. 
a distances of ~3-5 primary radii) in nearly circular or- 
bits and co-planar with the (new) rotational axis of the 
primary. For Haumea, the formation of the satellites is 
presumably concurrent with the formation of the fam- 
ily billions of years ago (jRagozzine fc Brownl l200l and 
the satellites have undergone significant tidal evolution 
to reach their current orbits (B05). 

4.3.1. Tidal Evolution of Semi-major Axes 

The equation for the typical semi-major axis tidal 
expansion of a single-sat ellite due to primary tides is 
(|Murrav fc Dermotti[20M) : 
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where ki p is the second-degree Love number of the pri- 
mar y, Q v is the primary tidal dissipation parameter (see, 
e.g.. iGoldreich fc Soterlll966l ). q is the mass ratio, R p is 
the primary radius, a is the satellite semi-major axis, and 

n = J GA ^ ot is the satellite mean motion. As pointed out 

by B05, applying this equation to Hi'iaka's orbit (using 
the new-found q = 0.0045) indicates that Haumea must 
be extremely dissipative: Q p ~ 17, averaged over the age 
of the solar system, more dissipative than any known ob- 
ject except ocean tides on the present-day Earth. This 
high dissipation assumes an unrealistically high k 2 ~ 1.5, 
which would be achieved only if Haumea were pe rfectly 
fluid. Using the strength of an rocky body and the lYoden 
(1995) method of estimating k 2 , Haumea's estimated k 2p 
is ~ 0.003. Such a value of k 2p would imply a absurdly 
low and physically implausible Q p <C 1. Starting Hi'iaka 
on more distant orbits, e.g. the current orbit of Na- 
maka, does not help much in this regard since the tidal 
expansion at large semi-major axes is the slowest part 
of the tidal evolution. There are three considerations, 
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however, that may mitigate the apparent requirement of 
an astonishingly dissipative Haumea. First, if tidal forc- 
in g creates a tid al bulge that lags by a constant time (as 
in Mignard 1980), then Haumea's rapid rotation (which 
hardly changes throughout tidal evolution) would natu- 
rally lead to a significant increase in tidal evolution. In 
other words, if Q p is fr equency dependent (as it se ems to 
be for solid bodies, see lEfroimskv fc L ainev 2007)), then 
an effective Q of ^16 may be equivalent to an object 
with a one-day rotation period maintaining an effective 
Q of 100, the typically assumed value for icy solid bodies. 
That is, Haumea's higher-than-expected dissipation may 
be related to its fast rotation. Second, the above calcu- 
lation used the volumetric radius R p ~ 650 km, in cal- 
culating the magnitude of the tidal bulge torque, where 
we note that Equation [2] assumes a spherical primary. A 
complete calculation of the actual torque caused by tidal 
bulges on a highly non-spherical body is beyond the scope 
of this paper, but it seems reasonable that since tidal 
bulges are highly distance-dependent, the volumetric ra- 
dius may lead to an underestimate in the tidal torque and 
resulting orbital expansion. Using R p ~ 1000 km, likely 
an overestimate, allows k2 P /Q p to go down to 7.5 x 10~ 6 , 
consistent both with k 2p — 0.003 and Q p ~ 400. Clearly, 
a reevaluation of tidal torques and the resulting orbital 
change for satellites around non-spherical primaries is 
warranted before making assumptions about the tidal 
properties of Haumea. 

The third issue that affects Hi'iaka's tidal evolution 
is Namaka. While generally the tidal evolution of the 
two satellites is independent, if the satellites form a res- 
onance it might be possible to boost the orbital expan- 
sion of Hi'iaka via angular momentum transfer with the 
more tidally affected Namaka. (Note that applying the 
semi-major axis evolution questions to Namaka requires 
a somewhat less dissipative primary, i.e. Q p values ^8 
times larger than discussed above.) Even outside of reso- 
nance, forced eccentricities can lead to higher dissipation 
in both satellites, somewhat increasing the orbital expan- 
sion rate. Ignoring satellite interactions and applying 
Equation [2] to each satellite results in the expected rela- 
tionship between the mass r atio and the semi-major a xis 
ratio of (jCanup et alj|1999t iMurrav fc Dermottll2000D : 

/ \ 13/2 
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where evaluating the right-hand side using the deter- 
mined orbits implies that — ~ 75.4 ±0.4. This mass ra- 
tio is highly inconsistent with a brightness ratio of ~3.7 
for the satellites, implying that the satellites have not 
reached the asymptotic tidal end-state. Equation [5] is 
also diagnostic of whether the tidally-evolving satellites 
are on converging or diverging orbits: that the left-hand 
side of Equation [3] is greater than the right-hand side im- 
plies that the satellites are on convergent orbits, i.e. the 
ratio of semi-major axes is increasing and the ratio of the 
orbital periods is decreasing (when not in resonance). 

4.3.2. Tidal Evolution of Eccentricities and Inclinations 

Turning now from the semi-major axis evolution, we 
consider the unexpectedly large eccentricities and non- 
zero inclination of the Namaka and Hi'iaka. None of the 
aforementioned considerations can explain the highly ex- 
cited state of the Haumea satellite system. As pointed 



out by B05, a simple tidal evolution model would require 
that the eccentricities and, to some extent, inclinations 
are significantly damped when the satellites are closer to 
Haumea, as eccentricity damping is more efficient than 
semi-major axis growth. One would also expect that the 
satellites formed from a collision disk with low relative 
inclinations. So, while a mutual inclination of 13° is un- 
likely to occur by random capture, a successful model for 
the origin of the satellites must explain why the satellites 
are relatively far from co-planar. 

The unique current orbital state of the Haumea sys- 
tem is almost certainly due to a unique brand of tidal 
evolution. Terrestrial bodies are highly dissipative, but 
none (except perhaps some asteroid or KBO triple sys- 
tems) have large significantly interacting satellites. On 
the other hand, gas giant satellite systems have multi- 
ple large interacting satellites, but very low dissipation 
and hence slow semi-major axis change. The change 
in semi-major axes is important because it causes a 
large change in the period ratio, allowing the system to 
cross many resonances, which can strongly change the 
nature of the system and its evolution. Even though 
Haumea is in a distinct niche of tidal parameter space, 
we can gain insights from studies of other systems, such 
as th e evolution of satellites in the Uranian system 
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or the interactions of 
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tidally-evolving exoplan ets (e.g.. IWu fc Goldreichl 120021 
iFerraz-Mello et ai1l2003l ). In addition, since the results 
of some Moo n-forming impac t s resulted in the cr eation 
of two moons (|Ida et al.lll997t) . lCanup et all (|1999D stud- 
ied the tidal evolution of an Earth-Moon-moon system, 
in many ways similar to the Haumea system. 

Using the results of these former investigations, we can 
qualitatively explain the excited state of the Haumean 
system. As the satellites were evolving outward at dif- 
ferent rates, the ratio of orbital periods would periodi- 
cally reach a resonant ratio. For example, as the current 
system is in/near the 8:3 resonance, it probably passed 
through the powerful 3:1 resonance in the relatively re- 
cent past. Since the satellites are on convergent orbits, 
they would generally get caught into these past reso- 
nances, even if their early eccentricities and inclinations 
were low. Note that this simple picture of resonance 
capture must be investigat ed numerically (as in, e.g., 
iTittemore fe Wisdom! 119881) and that higher-order res- 
onances may act dif f erentl y from lower-order resonances 
(| Zhang fc Hamilton! I2008D . Further semi-major axis 
growth while trapped in the resonance rapidly pumps 
eccentricities and / or inclinations, depend ing on the type 
of resonance (e.g.. IWard fc Canupl 12006). This contin- 
ues until the satellites chaotically escape from the reso- 
nance. Escaping could be a result of either excitation (re- 
lated to secondary resonances, see ITittemore fc Wisdom! 
Il989t iMalhotra fc Dermottlll990f) and/or chaotic insta- 
bility due to overlap ping sub-resonances, split by the 
large J 2 of Haumea ( Dermott et al.|[l988h . Outside of 
resonances, tidal dissipation in the satellites can damp 
eccentricities while the satellites are close to Haumea, 
but at their current positions, eccentricity damping is 
very ineffective even for highly dissipative satellites. In- 
clination damping is generally slower than eccentricity 
damping and was probably small even when the satel- 
lites were much closer to Haumea. Therefore a "recent" 
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excitation by passage through a resonance (possibly the 
3:1) can qualitatively explain the current orbital config- 
uration, which has not had the time to tidally damp to 
a more circular co-planar state. Numerical integrations 
will be needed to truly probe the history of this intrigu- 
ing system and may be able to constrain tidal parameters 
of Haumea, Hi'iaka, and/or Namaka. 

Note that early in the history of this system when the 
satellites were orbiting at much smaller semi-major axes 
(a < 10-Rp), the tri-axial nature of the primary would 
have been much more important and could have sig - 
nificantly affected the satellite orbits (|Scheeresl 11994) . 
At some point during semi-major axis expansion secular 
resonances (such as the evection resonance) could also 
be im portant for exciting ecce ntricities and/or inclina- 
tions (|Touma fc Wisdomlll998l ). Finally, tides raised on 
Haumea work against eccentricity damping; for certain 
combinations of primary and satellite values of A^/Q, 
tides on Haumea can pump eccentricity faster than it is 
damped by the satellites (especially if Haumea is partic- 
ularly dissipative) . While eccentricity-pumping tides on 
Haumea may help in explaining the high eccentricities, 
producing the large mutual inclination (which is hardly 
affected by any tidal torques) is more likely to occur in a 
reso nance passage, as with Miranda in the Uranian sy s- 
tem (jTittemore fc Wisdom|[T989HDermott et al.lll988h . 

5. CONCLUSIONS 

Using new observations from the Hubble Space Tele- 
scope, we have solved for the orbits and masses of the 
two dynamically interacting satellites of the dwarf planet 
Haumea. A three-body model, using the parameters with 
errors given in Tabled provides an excellent match to the 
precise relative astrometry given in Table [TJ The orbital 
parameters of Hi'iaka, the outer, bri ghter satellite matc h 
well the orbit previously derived in iBrown et al.l (|2005f ) . 
The newly derived orbit of Namaka, the inner, fainter 
satellite has a surprisingly large eccentricity (0.249 ± 
0.015) and mutual inclination to Hi'iaka of (13.41° ± 
0.08°). The eccentricities and inclinations are not due to 
mutual perturbations, but can be qualitatively explained 
by tidal evolution of the satellites through mean-motion 
resonances. The precession effect of the non-spherical 
nature of the elongated primary, characterized by J2, 
cannot be distinguished from the precession caused by 
the outer satellite in the current data. 

The orbital structure of Haumea's satellites is unlikely 
to be produced in a capture -related formation mechanism 
(e.g.. IGoldrcic rTet al.ll2002f L Only the collisional forma- 
tion hypothesis (allowing for reasonable, but atypical, 
tidal evolution) can explain the nearly co-planar satellites 
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(probably with low den sities), the rapid spin and elon- 
gated shape of Haumea (Rabino witz et ai1l2006l ). and the 
Haumea collision al family of icy obj ects with similar sur- 
faces and orbits (|Brown et al 1 [20071) . 

The future holds great promise for learning more about 
the Haumea system, as the orbital solution indicates that 
Namaka and Haumea are undergoing mutual events for 
the next several years. This will provide excellent obser- 
vational constraints on the size, shape, spin pole, density, 
and internal structure of Haumea and direct measure- 
ments of satellite radii, densities, and albedos. There 
are also interesting avenues for future theoretical inves- 
tigations, especially into the unique nature of tidal evo- 
lution in the Haumean system. These insights into the 
formation and evolution of the Haumean system can be 
combined with our understanding of other Kuiper belt 
binaries to investigate how these binaries form and to 
further decipher the history of the outer solar system. 
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